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Abstract 

The p^p^ decay mode is searched for in a data sample of about 227 million T(45) 

BB decays collected with the BABAR detector at the PEP-II asymmetric B factory at SLAC. No 
significant signal is observed, and an upper limit of 1.1 x 10~^ (90% C.L.) on the branching fraction 
is set. Implications on the penguin contribution and constraints on the CKM angle a with B ^ pp 
decays are discussed. All results are preliminary. 
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1 INTRODUCTION 



Measurements of CP-violating asymmetries in the B^B^ system allow tests of the Standard Model 
by over-constraining the Unitarity Triangle through the measurement of its angles. The time- 
dependent CP asymmetry in a 6 — > uud decay of a B^ to a CP eigenstate allows for a direct 
measurement of the angle a = arg [—ytd^tb/^ud^ub\ decay is dominated by the tree amplitude. 
The contribution from penguin diagrams gives rise to a correction Aq = Ucq — ol that can be inferred 
through an isospin analysis as illustrated in Fig. Q The theoretical error from isospin-breaking 
effects is typically estimated to be smaller than 10 degrees |2]. 




I = 1 

Figure 1: Isospin Triangles relating the amplitudes of the modes B^ P -, ~^ ■, ~^ 

p~^p^ and their charge conjugates for a given polarization (e.g., CP-even longitudinal polarization). 

The 7r"'"7r~ final state (CP =1) was expected to be an ideal candidate for the study of a. However, 
the recent measurement of B^ — > vr'^vr'^ confirms that this mode suffers from substantial penguin 
contamination [2]. This severely limits the possibility of a model-independent measurement of a 
with currently available statistics. On the other hand, the recent experimental results shown in 
Table ^ confirm the theoretical expectation |^ of a smaller penguin contribution in the pp system 
with respect to the vrvr case (charge-conjugate decay modes are always assumed throughout this 
paper). 

Table 1: Recent measurements of B ^ pp branching fractions. 



Decay BABAR Belle 



BR{B+ - 


- P^P': 


) X 10^ 


22.5l|^±5.8 31.7±7.ltlj M 


BR{B° - 


P'^P" '. 


) X 10^ 


30±4 ±50 


BR{B^ - 




X 10^ 


<2.1 (90 % C.L.) 



In pp decays the final state is not, in general, a CP eigenstate, and an isospin triangular 
relation holds for each of the three helicity states, which can be separated through an angular 
analysis. However, the measured polarizations in p'^p~ and p~^p^ modes indicate a dominance of 
the helicity state (longitudinal polarization), that is, of a CP =+1 eigenstate. A measurement 
of the polarization in B^ p^p^ would complete the isospin triangle, but this mode has not been 
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observed so far. The best present limit on the decay was obtained by BABAR with a sample 
of 89 million T(4S') BE decays p]- The resulting Grossman-Quinn |H] bound that can be set 
on the penguin contribution [71 IH] is more stringent than in the vrvr case, but knowledge of the 
B^ — > p^ p^ rate is still expected to be the limiting factor to the accuracy of the a measurement 
with pp decays. In jS], Grossman and Quinn show that the angle 2Aa^ (see Fig. ^ between the 
two triangles is maximum when the two isospin triangles have opposite direction and they are right 
triangles. In that case, the bound takes simply the form 

where and f^' are the fraction of longitudinal polarization in B'^ p^p^ and B^ — > p^p*^ 
decays respectively, and B{B^ p'p') and B{B~^ P^p') the respective branching fractions. 

In this paper, we present a search for the p^p^ final state performed at BABAR on a sample of 
227 million T(AS) — > BB decays. Improvements over the previous result are achieved from the 
increased statistics and from optimizations of the analysis, which result in an increased sensitivity. 



2 THE BABAR DETECTOR AND DATASET 

The BABAR detector operates at the PEP-II asymmetric B Factory at SLAG and is described in 
detail elsewhere jlUj . The data samples used in this analysis correspond to an integrated luminosity 
of 205.4 fb"^ collected at the T{AS) resonance, and of 16.1 fb"^ cohected about 40 MeV below 
the T{4:S) resonance in order to study the qq continuum background. Large samples of Monte 
Garlo data are also used to model signal B backgrounds. 



3 ANALYSIS METHOD 

We fully reconstruct B^ p^p^ candidates from their decay products p^ vr'^vr" with four charged 
tracks in the final state. Gharged track candidates are required to originate from a single vertex 
near the interaction point. We select p*^ candidates with requirements on the ■k~^tt~ invariant mass, 
loose enough to retain sidebands for later fitting. 

Improvements in the track reconstruction efficiency and in the understanding of the associated 
systematic error have allowed us to extend the accepted momentum range with respect to the 
previous analysis, and consequently increase the selection efficiency by 25%. The particle identi- 
fication capabilities of the BABAR detector are used to reject tracks identified as electrons, kaons, 
or protons. A set of kinematic variables describing the shape of the event is used to suppress the 
qq continuum background. The identification of signal events is based on two kinematic variables: 
the beam-energy substituted mass of the B 

rriES = \l{s/2 + • PbY/E^ - p| (2) 

(where the initial four-momentum (£^j,pi) and the B momentum pB are defined in the laboratory 
frame), and the difference between the reconstructed B energy in the center-of-mass frame and its 
known value 

= E'i^ - (3) 
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An unbinned extended maximum likelihood fit is then performed on the selected sample to evaluate 
the signal yield using several discriminating variables. The probability density functions (PDF) 
entering the likelihood function are built as the sum of several components that describe separately 
the signal, the qq continuum, and various background BB decays. The result was revealed only 
after the finalization and validation of the full analysis. 

3.1 Angular Observables 

The angular distribution of the B meson decay to a vector-vector iVV) final state is a combination 
of S-, P-, and D-wave contributions with unknown relative amplitudes. The helicity angles (^j, 
(/) = (^1 — (/>2, i = 1, 2) are defined by the direction of the two-body decay axis and the direction 
opposite the B in the rest system, as shown in Fig. |21 




Figure 2: Definition of helicity angles ^i, 62, and 0, for the decay B — > pp. The p final states 
are shown in their rest frames. 



Due to quark spin conservation, the longitudinal polarization is expected to dominate the 
VV decays. However, this expectation might be strongly violated in penguin-dominated decays, as 
has been observed m. B ^ (j)K* (Hj. Since the azimuthal angle (p does not provide any significant 
background suppression, and there are no acceptance effects for this observable, we integrate over 
it, which results in: 



1 



r d cos 9i d cos 62 



9 

8^ 



cos 



h cos^ + 7 sin2 01 sin2 62 (1 - fi) 



(4) 



The two helicity observables 7ii = \ cos 9i \ are used in the fit to measure the longitudinal 
polarization of the decay, fi = |Aop/(|Aop + |^+ip + |A_ip), where A\ (A = 0, +1,— 1) correspond 
to the helicity amplitudes of the two /o's in the B rest system. These observables also provide 
additional background suppression. The limit on the branching fraction can be expressed as a 
function of fL- The hypothesis fL = ^ gives the most conservative upper limit on the amplitude, 
shown in Fig. ^ This is also the relevant result for the measurement of a through the isospin 
analysis. 
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3.2 Event Selection 



In order to reject the dominant quark-antiquark continuum background, we require | cos^tI < 0.8, 
where Ot is the angle between the i?-candidate thrust axis and that of the remaining tracks and 
neutral clusters in the event, calculated in the CM frame. The other event-shape discriminating 
variables include the polar angles of the B momentum vector and the S-candidate thrust axis with 
respect to the beam axis in the T(4S') frame, and the two Legendre moments Lq and L2 of the 
energy flow around the i?-candidate thrust axis. These variables are combined in a neural network, 
which is trained using off-resonance beam data and signal Monte Carlo data. The neural network 
output is transformed into a Gaussian-shaped variable (called iS-variable hereafter), without any 
loss of discriminant power, to make the shape easier to fit. 

To further suppress background, we use multivariate algorithms to identify the flavor of the 
other B in the event (tagging) The suppression power comes from the fact that signal and 

background have different tagging efficiencies in the various tagging categories. These categories 
correspond to different methods of identifying the h content of the other B meson in the event. We 
use five tagging categories Ctag (lepton, two kaon, inclusive, and the untagged event categories). 
Since BB events tagged by the best performing tagging categories (notably by leptons) have a lower 
qq background, the splitting into tagging categories improves the discrimination of the continuum. 

To summarize, we use eight primary observables to characterize each event candidate: Xi = 
{mEs, AE', <5, m(7r+7r~)i, m(7r"'"7r~)2, TYi, 'H21 Ctag}- The ranges for the key discriminating variables 
are the following: we require tties > 5.24 Gev/c^ and |Ai?| < 85 MeV (the signal mss ^-nd 
resolutions are 2.5 MeV/c^ and 20 MeV respectively), m(7r"*'7r~)j between 0.55 and 1.0 GeV/c^, 
and 7ii < 0.99. The latter cut removes a region with very low efficiency. 

Additional selection criteria are finally applied to veto the most dangerous source of BB to 
charm background, namely B^ D~Ti^ — > if"'"(7r"*")7r~7r~7r+, by requiring the invariant mass of 
the three-particle combination that excludes the highest-momentum track in the B frame to be 
inconsistent with a D meson. 

After applying all selection criteria, the event sample has on average 1.05 candidates per event. 
When several candidates are selected in the same event, one candidate is selected randomly. The 
selection efficiency, measured on signal Monte Carlo samples, is 27% (32 %) for the longitudinaly 
(transversely) polarized signal. This corresponds to an expected signal yield in our data sample 
of about 60 events for a branching fraction of 1 x 10~^. The fraction of selected events where 
the four candidate tracks do not match the true event tracks, called self-cross-feed (SCF) events 
in the following, is 22% (8%) for longitudinal (transverse) polarization. The large self-cross-feed 
for longitudinal signal is due to the presence of soft pions that can be easily exchanged with soft 
particles from the other B. 
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3.3 Backgrounds 

The selected sample is expected to be dominated by qq combinatorial background. However, the 
most dangerous backgrounds are expected to come from other B decays, in particular from final 
states with four charged particles that can mimic the signal in the distribution of uies and AE, 
for four-pion final states. As mentioned above, the B to charm decays are suppressed by explicit 
veto cuts. For the B to charmless decays, we expect contaminations from more than 50 modes, the 
most important ones being listed in Table [21 The mode B^ — > afir^ — > p'^vr^vr^ is an irreducible 
source of background. The discrimination of this mode relies mainly on the p mass distributions. 

The other charmless modes have smaller impact on the result. The most important ones are 
described by specific PDFs, in order to improve the background model in the likelihood fit. We 
fix the yields of well measured modes {p~^p^, p'^P~, P"^) and modes which are isospin-related to 
measured modes (p^K*^), and vary these yields to evaluate the corresponding systematic error. 



Table 2: Background categories which are either floated or fixed in the fit to the data sample. The 
yield uncertainties in the fixed background categories (both statistical and systematic) are taken 
into account in the evaluation of the systematic error. 



Background Category 



Yield 



B - 

B^ 

B^ 

B+ 

B - 

B° 



charm 

-.p+p' 
p-K modes 



25 



other charmless B decays 



floated 
floated 
floated 
± 5 ± 16 
78 ± 9 ± 18 
34 ± 6 ± 4 
14 ± 4 ± 3 
floated 



m 
m 



3.4 Likelihood Fit 

The signal is obtained by maximizing the extended likelihood function 

(A^cat \ A^cand / \ 
-E™^- I n (E n,V,{xfJ)\ (5) 

with several event categories j: signal (including a SCF fraction), continuum qq, and several BB 
background categories listed in Table |2 Each Vj{xi;(3) is the PDF for the observable Xj, and is 
described by the PDF parameters (5. 

The event yields rij in different categories are obtained by minimizing the quantity — 21n£. 
The statistical significance of a signal is defined as the square root of the change in — 21n£ when 
constraining the number of signal events to zero in the likelihood fit. The assumption of negligible 
correlations among most of the fit input variables has been validated with MC simulation and data. 
However, there are effects which do introduce noticeable correlations in the samples of signal and 
background events, such as helicity angles correlation in signal, and mass-helicity correlation in 
background, as discussed below. 
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3.5 PDF Parameterization 

We use double-Gaussian functions to parameterize the AE and mss PDF's for signal, and a 
relativistic P-wave Breit-Wigner convoluted with a Gaussian resolution function for the resonance 
masses. The £^-variable is described by an asymmetric Gaussian plus a single Gaussian function for 
both signal and background. The signal helicity PDF is expressed as a function of the longitudinal 
polarization (see Eq. The ideal angular distribution is multiplied by the detector acceptance 
function G{'Hi,'H2)- We obtain the acceptance function from a fit to a sample of MC events with 
transverse and longitudinal polarization. 

For the combinatorial background we use low-degree polynomials for AE and resonance masses, 
and an empirical phase-space function for rriES (known as Argus function jl3j): 

f{x) oc — ^2 exp[— ^(1 — x^)] , (6) 

where x = mEs/Ehcam, and ^ is a parameter measuring the curvature of the distribution near the 
end-point. The background parameterization of resonance masses also includes a resonant com- 
ponent to account for resonance production, which is described by the same P-wave Breit-Wigner 
function used for the signal p shape. The background helicity-angle distribution is separated into 
contributions from combinatorial background and from resonances. The resonances are assumed 
to be unpolarized. We parameterize the combinatorial helicity distribution with a second-degree 
polynomial and an exponential function to allow for the increased fraction of fake p candidates 
with low momentum pions near TCi = 1. The amount of peaking near Tii = 1 depends on the p 
candidate mass and is parameterized accordingly. Fig. |31 shows the m-Es, AE, and £ PDFs for 
combinatorial background and longitudinal signal. The two-dimensional mass-helicity distribution 
for combinatorial background, and the projection of the helicity distribution for longitudinal signal 
are also shown. 

The PDFs for exclusive B decay modes are modeled with non-parametric distributions describ- 
ing the shapes of the observables, except when certain distributions are expected to be identical to 
those used for signal. These distributions are described by smoothed histograms ^1] with a large 
number of bins. In the S-background modes the two p^ candidates can have very different mass 
and helicity distributions. This happens, for example, when one of the two p^s is real (e.g. p'^p^, 
p^K*^) or when only one of the two p^s contains a hard bachelor pion (oivr). In such cases, we con- 
sider a four- variable correlated mass/helicity PDF. The inclusive i?-to-charm and i?-to-charmless 
categories are parameterized similarly to the exclusive modes, they include the remaining B decay 
modes not modeled explicitely, and their yield is left free in the fit, to include possible unaccounted 
background sources and compensate for possible imperfections in the background model. 

The B-Havor tagging PDFs for signal and background are simply the discrete distributions of 
tagging efficiencies. Large samples of fully reconstructed B meson decays are used to obtain the 
-B-tagging efficiencies for signal B decays and to control the MC values of S-tagging efficiencies for 
the B backgrounds. Continuum background efficiencies are obtained from the sideband data. 
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Figure 3: PDFs of the observables (from top to bottom): rriEs? Ai?, and for combinatorial back- 
ground (left) and longitudinal signal (right). The bottom row shows the combinatorial background 
mass-helicity two-dimensional PDF (left), and the projection of the helicity distribution for signal 
(right). The small dip near | cos^| = 0.85 is due to the D-veto acceptance effect. 
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4 PHYSICS RESULTS 



Table |31 shows the results of the fit. No significant yield is observed, and an upper limit on the 
branching fraction of — > p^p'^ is set. As mentioned earlier, 100% longitudinally polarized signal 
is used in the fit, as it was checked that it gives the most conservative estimate for the upper limit. 

Table 3: Results of the fit: Signal yield {Ns), selection efficiency (Eff), branching fraction 
(i3), upper limit at 90% confidence level (U.L.) and significance of the measurement, defined as 
y^A(— 21nvC) when constraining the number of signal events to zero in the fit. The first error cor- 
responds to the statistical uncertainty and the second one to the systematic uncertainty, discussed 
in next section. 

Quantity Value 

Ns 

EflF (%) 

i3(xl0-6) 
U.L.(xlO"6) 

Significance (y^A(— 21n£)) 



Figure m shows the result of the fit projected onto the thes s^nd AE observables. The histograms 
show the data after a cut on the quantity Vsig/ ("Psig + 'Pback) has been applied, where Vsig and T-'back 
are the probabilities for a given event to be signal and background respectively, and are evaluated 
using all the observables except the one that is being plotted. The cut is optimized for each variable 
separately. The solid (dashed) line shows the projection for the full fit (background only) after the 
same cut is applied. 



331^5 ± 12 
27.1 ± 1.3 
0.54l°J^ ±0.19 
1.1 (1.0 statistical only) 
1.6 (1.9 statistical only) 
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Figure 4: Fit result projections onto niEs and IS.E. The histograms correspond to the data and the 
sohd (dashed) Hne to the full (background only) fit after a cut on the quantity 'Psig/('Psig + "Pback) 
is applied. The projections contain 22.5% and 23.9% of signal, and less than 0.5% and 0.2% of 
continuum background respectively. 
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5 SYSTEMATIC STUDIES 



Table |^ summarizes the systematic errors on the branching fraction. We include the uncertainty 
in the total number of B mesons, the uncertainties coming from selection cuts, such as track 
multiplicity, thrust angle, and vertex requirements, and the uncertainty in particle identification. 
The uncertainty on the selection efficiency is estimated from Monte Carlo and is dominated by the 
reconstruction of soft tracks. The accuracy of the simulation of the key observables mss and AE 
is estimated using the — > D~-k~^ — > tt'^ control sample. 

The peaking BB background has been taken into account in the fit with three fixed components. 
The effect of their systematic uncertainty on the signal is evaluated by adding to the data known 
samples of background Monte Carlo events and observing the variation of the result. 

Monte Carlo studies show that the interference between signal and a^vr''^ may not be neglected 
and could significantly bias the measurement. Assuming for the branching fraction of the two modes 
the central values measured without accounting for interference, and accounting for acceptance and 
likelihood fit, the systematic bias is estimated at the level of 7.5 events. 

To obtain the error associated to the uncertainty in the PDFs parameters, 200 Monte Carlo 
experiments are performed with the parameters varied within their errors. The width of the fitted 
yield distribution is taken as systematic error. Finally, we assign a systematic error of it 3 events 
to cover a possible fit bias evaluated with Monte Carlo experiments. 

Table 4: Summary of systematic errors in the measurement of the B^ p^p^ branching fraction. 



Source 


Uncertainty 




Yield fraction Number of events 


Multiplicative 



1.1% 0.4 

1% 0.3 

1% 0.3 

2% 0.7 

2% 0.7 

3.2% 1.1 

<1% <0.3 



Additive 


B background 


17.4% 


5.8 


aivr interference 


22.5% 


7.5 


PDF variation 


18% 


6.0 


Fit bias 


9% 


3.0 


Total 


35% 


11.7 



Taking systematic uncertainties into account, the measured value for the B p p branching 
fraction is 

B{B p^p^) = (0.54;°:| ± 0.19) X 10^*5 
or an upper limit of 1.1 x 10~^ at 90% confidence level. 



Number of B mesons 
Track multiplicity cut 
Thrust angle cut 
Vertex requirement 
PID cut 
Track finding 
MC statistics 



17 



6 SUMMARY 



With a sample of 227 million T{4S) — > BB decays collected with the BABMi detector we have 
searched for the decay mode B^ p^p^. We set an upper limit of 1.1 x 10~^ (90% C.L.) on the 
branching fraction of this decay mode. Our results are preliminary. This result has important 
implications for understanding the penguin contribution and constraints on the CKM angle a with 
the B ^ pp decays. With the B^ P^P~ values of the branching fraction and longitudinal polar- 
ization [7j, the measured B^ — > p'^p~ S and C CP-violating time-dependent asymmetry parame- 
ters jni; the B^ P^P^ values of the branching fraction and longitudinal polarization [2] [S] |12j . 
and the measured branching ratio of B^ p^p^ measured in this analysis, neglecting electroweak 
penguins, non-resonant and 1 = 1 isospin contributions, and using isopin analysis |15| . we obtain a 
new value for a: 

a = (96 ± lO(stat) ± 4(syst) ± ll(penguin))° 

To extract alpha a i^i which all the measured quantities and the angles 2Aa^ and 2Aa^^ 
of Fig. Hare expressed as a function of the length of the two isospin triangles' sides is mimimized. 
The confidence level on a is obtained by a scan of the difference between x^i^i) for a given value 
of Q and the minimum of this x^- Fig. |S1 presents the result of this scan. 




a (deg) 



Figure 5: Confidence Level on a obtained from the measured branching fraction, fraction of lon- 
gitudinal polarization, and CP parameters of B^ — > p~^p~, branching fraction and fraction of 
longitudinal polarization of B~^ — > p'^p^, and the B^ p^p^ branching fraction measurement 
from this analysis, using isospin analysis. The dashed lines correspond to the 68% (top) and 90% 
(bottom) confidence intervals. 
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